Channel coupling in heavy quarkonia: energy levels, mixing, widths and new states 
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The mechanism of channel coupling via decay products is used to study energy shifts, level 
mixing as well as the possibility of new near-threshold resonances in cc, bh systems. The Weinberg 
eigenvalue method is formulated in the multichannel problems, which allows to describe coupled- 
channel resonances and wave functions in a unitary way, and to predict new states due to channel 
coupling. Realistic wave functions for all single-channel states and decay matrix elements computed 
earlier are exploited, and no new fitting parameters are involved. Examples of level shifts, widths 
and mixings are presented; the dynamical origin of X(3872) and the destiny of the single-channel 
2 3 Pi(cc) state are clarified. As a result a sharp and narrow peak in the state with quantum numbers 
, J PC — 1 ++ is found at 3.872 GeV, while the single-channel resonance originally around 3.940 GeV, 

becomes increasingly broad and disappears with growing coupling to open channels. 

PACS numbers: 12.39.-x,13.20.Gd,13.25.Gv,14.40.Gx 
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Most hadron states are coupled by strong interaction to closed or open decay channels, and thus are subjects of 
the Theory of Strongly Coupled Channels (TSCC). The latter topic was developed during many decades, see [Tj-ll] 
and [5] for a review, and also as more recent publications. In the present paper we apply TSCC specifically to 
the case of Okubo-Zweig-Iizuka rule allowed two-body decay channels of charmonia and bottomonia. In doing so we 

^ . need several prerequisites. First of all ^ it is the one-channel description of charmonia and bottomonia as cc, bb states 

in relativistic Hamiltonian formalism Q , developed in the framework of the Field Correlator method (FCM) [13] (see 
flTj [ for a review) and wave functions of stable states in x or p space cast in the numerical form. The latter have 
been accurately computed using this method with only universal input: the string tension a, the current (pole) quark 
masses m,i, and the strong coupling a s (q) [l2j . 
£ — ■ The next ingredient is an effective relativistic Lagrangian for the pair creation, inducing the string breaking. To 
this end we are using the decay mass vertex J -ipM^tp d 4 x introduced in [l3j and exploited for dipion transitions in 
[l3T - [l5j and for the reaction channel T(nS) — > BB, BBn in [ll|. In principle, M w can be expressed in terms of quark 
masses and average energies, but we use it as the only one parameter, which is fixed in our previous studies [l3l - fl6j 
J> , Finally, as it was shown in [l3[ and before in [13, [H| the transition matrix element reduces to the overlap integral of 
wave functions of decaying system and products of decay. It is interesting, that the vertex operator in this integral 
contains not only M w , but also the Z factor of the decay process constructed from the Dirac trace of all involved 
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hadron vertex states, and projection operators. This technic, introduced in [131 ] ■ is a relativistic equivalent of the 
nonrelativistic one with spin-angular momentum (Clebsch- Gordon) coefficients used in the framework of 3 Pq model 
flfjj . As a result one obtains a system of integrodifferential equations for new wave functions and energy eigenvalues, 
which can be easily solved in the lowest approximation for energy shifts, widths and level mixing coefficients. At the 
same time we have developed a (2 x 2) variant of wave functions and matrix elements for light quarks in heavy-light 
mesons. Several examples of this kind are shown below. 

At a deeper level, one meets with several problems: (i) First, the states above decay thresholds are unstable and the 
definition of the wave function itself is questionable in a rigorous sense, since an admixture of continuous spectrum 
states appears. Here different approaches exist. The most rigorous is the Weinberg procedure [201 ] . named below the 
Weinberg Eigenvalue Method (WEM). It is used to define the resonance wave function and energy, as well as t- matrix 
via Weinberg eigenvalues. The WEM has been used before in the one-channel 2-body and 3-body problems [20l| . [2lj| . 
We have found that it is specifically useful in case of TSCC, since coupled channels (CC) induce the energy-dependent 
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force term, which violates standard orthonormalization procedure, while this term can easily be treated in WEM. (ii) 
Second, even the closed channels cause the problems. In terms of hadron loops it was treated in many papers, see e.g. 
and recent papers 0, H[ , where some theorems were formulated 0] and the renormalization method was suggested 

In essence the problem here is similar to the problem of unquenched quark pairs. It occurs also for stable hadrons, 
where the renormalozation procedure is necessary in general. We will not discuss these topics in the given paper, 
assuming that the renormalization is done e.g., by readjusting the pole quark mass. We also disregard the important 
topic of full relativistic invariance for composite objects moving with different velocities, e.g. charmonium decaying 
in its cm. system into heavy-light mesons, with their wave functions defined in their cm. systems. This is done 
assuming small relative velocities near thresholds. Far from thresholds these factors become important. Similarly, 
near thresholds we assume here, as well as in [l3l - [l6| . the 3 Pq type of the decay vertex, while at higher energies this 
type of decay may be replaced by another one, e.g. the 3 Si type. 

In this paper we systematically apply WEM to find the shifts and widths of (n 3 Si) energy levels, as well as mixing 
between them. We find the method to be especially useful to discover the analytic structure and pole positions in the 
case of strong CC. A particular example of the 2 3 Pi level proves to be a good illustration of our analysis. From the 
experimental point of view two interesting problems appear. First one, why in experiment only the peak at the lowest 
DqDq threshold is seen, while at the slightly higher, D + D*_ no peak was ever seen? Secondly, possible resonances 
at 3.940 GeV found in [22J, seemingly are not J = 1, and very likely the 1 ++ state around 3.940 GeV was never 
observed. Our analysis allows to answer both questions at the same time, as will be explained below. 

As a result we find two poles due to a single eigenvalue in the possitions near 3.872 GeV and 3.940 GeV, but the 
latter peak becomes too wide and finally disappears with increasing CC, which can explain the experimental situation 
[22j . Moreover, retaining the peak appears at the lower threshold 3.872 GeV, and not at the higher threshold 3.879 
GeV, again in agreement with experiment (22[. 

The plan of the paper is as follows. In the next section we introduce the general formalism for the Green's functions 
of charmonia and bottomonia with the inclusion of decay channels. We present equations for wave functions (Green's 
functions) both in QQ and {Qq){Qq) channels. In section Hill the CC resonances are considered in the decay channel 
and condition for the existence of a CC resonance is formulated. In section HVl the rigorous Weinberg theory of CC 
resonances is presented. In section [V] the mixing of states in WEM is considered. In section IVT1 we present results for 
values of level shifts and widths for 3 3 Si state and also mixing between 3 3 S\ and 2 3 Si state, as well as the analysis 
of the situation in the 2 3 Pi state. In section Ivnl summary and prospectives are given. 



II. GENERAL FORMALISM OF STRING-BREAKING CHANNEL COUPLING 



We consider two sectors of hidden and open flavor with initial and final bare gauge invariant operators, for heavy 
quarkonium sector, 

L 3i (I \ x ) = ^q{ x )^i ^q{x) 
and for heavy- light meson sector, 

ii. ti n \x) = $ Q (x) r, i> q (x), 

where Ti = 1, 7^, D^a^, .... With the help of j/ 7 ' 77 ^ one generates bare mesons and as shown in [l3], one can 
project physical amplitudes 1 (Green's functions) with physical wave functions , ^Qq \ ^Q^g ^ or stationary 

states one can use Green's functions in energy representation, e.g. 

G^(l,2;£)=y QqU QQ U = ^—. (1) 
QQ y ' ' > ^ E n , -E Hq-E y ' 

Here superscript (0) of Green's function refers to the bare case, when sector II is switched off, and ^qq, E ni refer to 
the eigenfunctions and eigenvalues of the relativistic string Hamiltonian Hq Q , for charmonium those were calculated 
in |l2| and for bottomonium in (25[. In sector II the counterpart of (JXJ) consists of Green's function of the pair 
(Qq), (qQ). We neglect in the first approximation interaction of two color singlet mesons, and write the cm. Green's 



1 Note, that the procedure of hadron state projection is here fully equivalent to that used in lattice approach, see e.g. [24j . 



FIG. 1: Decay matrix vertices. 



function as 

At this point we must take into account a possible transition (decay) of states in sector I into states of sector II, 
which can be done in several ways. In literature it is common to assume one of several types of phenomenological 
decay Lagrangians, e.g. 3 Pq type [l!|, with vector confinement vertex or with scalar confinement vertex, studied 
in [l8[ . For bottomonium a relativistic string decay vertex of the following form was used in [l3l - [l6| ] 

C sd = J 4>q ip q d A x, M u = const (3) 

where M w was taken to be constant, M u « 0.8 GeV from decays of bottomonium into BB, BB* , ... [l5l [l6|. 

It is important, that we work in the cm. system and consider both wave functions and Hamiltonian obtained in the 
instantaneous hyperplane, when all time coordinates of all particles are the same 2 . Therefore the vertex C s d enters 
between instantaneous wave functions of QQ on one side and product of Qq, qQ on another side 

J\2Z = ^= J Vl23^ qq M u i>Qqip q Q dr. (4) 

At this point one should define exactly how spin, momentum and coordinate degrees of freedom enter in Q , which 
we denote by an additional factor j/123 and yet undefined phase space factor dr. One way is to exploit nonrelativistic- 
type decomposition which is used in 3 Pq calculations [19j. In our considerations we are using two different ways. 
Below we begin from the fully relativistic formalism of Dirac traces and projection operators, started in [23] for decay 
constants and in fl3l ] for dipion transitions and then will go into reduced (2 x 2) spin-tensor formalism, explained in 
Appendix B. Note that the relativistic formalism with Z factors is similar in lattice calculations for transition matrix 
elements, see e.g. [24j ], 

In the formalism one considers initial and final meson creation operators I, II, given above (see in Appendix B the 
Table IVT1 of lowest operators and their (2 x 2) forms), and composes the decay matrix element as shown in Figure 1 
where vertices 1,2, 3, a; are operators r i=1 .2,3.x entering in the bilinears ji — tpTitp, and lines (1, 2), (2, x), ... denote 
quark Green's functions 5q(1,2), Sq(2,x), etc. so that matrix element corresponding to Figure[T]is 

S{l,2,x,3)= tr S Q (l,2)T 2 S q (2,x)T x S q (x,3)T 3 SQ(3,l)Y (5) 

As was shown in [15, [231] in the approximation where one neglects influence of spin forces on wave functions, one can 
replace 



(mQ, q + ^Q,ql4 ~ ipf ,9 7-i) ( m Q,9 - w Q,g74 + ipf' q lt) 

S Q,9 = = A Q,q G Q,n S Qq = K— G Q,q = K Q,q G Q,q- ( 6 ) 



2 We omit boost corrections here, which makes application of our method justifiable only close to thresholds. Far from thresholds one 
should take into account both boost and, more important, a possible change in the pair creation vertex, since high energy transfer to 
the qq pair might require gluon exchange mechanism, hence 3 Si vertex 
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where Gq iQ is quadratic Green's function, Ag are projection operators, the variables WQ ig are the averaged kinetic 
energies and rag >q are the pole masses. 

Our physical matrix element corresponding to the decay — > ^Qq can obtained from ((SJ) by projecting 

the chosen intermediate states, as shown by dashed lines in Figure [TJ 

As a result as shown in [l3|, [l4j the physical projected matrix element has the form 

J m (p) = ^=J V%i *<$(u - v) eW^iu - x) 4"^(x - v) d 3 xd 3 (u - v) , (7) 

where N c is the number of colours, r = c (u — v), c = and j/123 = g • The expressions for Z123, 

are proportional to Dirac traces of the projector operators and are given in [13,14]. We point that the w.f 
^QQ 1 ^Q? ' ^Qg in ® are 110 longer full w.f. of mesons, but the radial part Rqq , Rq? , Rq* divided by 
while the angular part of the w.f. is accounted for in the factor j/ 12 3 (2/123)- Important role is played by average values 



of quark kinetic energies, wq,<j = (y m Q q + P 2 ) inside heavy-light mesons in their cm. systems (if one neglects cm. 

motion of these mesons). The numbers of ujq, uj q are computed from the relativistic string Hamiltonian in [T2I . [25[. 
In practical calculations it is more useful to exploit the (2 x 2) reduction of the bispinor wave functions in (|4]), as 

tjj = ( ^ j and ij> T ip = (v c , w c ) 7274 T ( ^ J , see details in Appendix B. The resulting matrix element has the same 
form as in ([7]) but with M w j/123 — > 72/123 1 where 2/123 is given in Table IVTll and 7 is proportional (see Appendix 

Now one can define the selfenergy part in sector I due to sector II in the intermediate state, which is 

/rnv f d S p v— * Jnn 2 n 3 (p) Jmn 2 n 3 (p) fa \ 

Wnm{E) = J W? £ ^Wp) (8) 

and the total Green's function in sector I can be written as (sum over bound states only) 
G W (12 . E) V *&(D*& W) (2) ^(1)^(^)^(2) 

where ellipsis implies terms of higher order in w nm and this can be summed up as 



n . n 1 



where matrix (E) nm = £"„ <5 nm . 

Note, that in (fTU| the Green's function is actually a projection of the coupled-channel system on the original 

unperturbed QQ wave functions ^qq- In reality wave functions of the coupled-channel system differ from the latter 
and acquire continuous spectrum pieces above the decay threshold, and hence need a special treatment to be discussed 
below. 

The new spectrum is obtained from (fTU| as 

det (E - E - w) = (11) 

for one level in sector I it simplifies 

E = E n +w nn (E) (12) 

which yields energy shift and width in the first order approximation in w: 

= E n + Re (w nn (E n )) , r« = 2 Im (w nn (E n )) . (13) 

In the next order one should solve the transcendental in E one-channel equation (|12[) . which is valid when w nn is 
large, but |io nm | < \E n - E m \. 

Below the decay threshold one can diagonalize the matrix in (|10[) with unitary matrices 

{{E-E- w)~\ m = [^-J— U Xm {E) (14) 
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and the Green's function acquires the form 

A n 

In this way § v become new orthogonal states comprising all effects of mixture between bound states due to closed 
channels. The same procedure can be applied for open channels (above the decay threshold) when one neglects the 
widths of the levels, i.e. imaginary part of w. 

One can define interaction V\2\ in sector I due to sector II, 

V 121 (r, r') = £ G^Jr - r') X mn3 (r) X+„ 3 (r') (16) 



with 



G (o) (r _ r ') = / ±_P_^Z_ (17) 

^nsW = ^=J |^ qr *WS?(<l) (18) 



or, in momentum space 



where X„ 2 „ 3 (Q) = ^= V'QgHQ) (Q)- Now the one-channel Hamiltonian i? in sector I is augmented by the 
term Viai, 

H = Hq + V 121 , HV ni = EV ni (20) 

Note, that the coupled-channel (CC) interaction can be strong enough to support its own bound states, as was studied 
in where this type of resonances was called the CC resonances. 

Let us now turn to the relativistic string Hamiltonian (RSH) Hq, derived from the gauge- invariant meson Green's 
function in QCD in the one-channel case in Q . This Hamiltonian has been successfully applied to light mesons (2(| , 
heavy- light mesons [23|, , and heavy quarkonia [IH, [25| and has a simple form: 

uji loo m? Too p 2 , . . . 

*° = T + f + 2^ + i + 2^ + ^ (21) 

V n (r)=V B (r)+V SD (r,uj i ). (22) 



In general, the quantity Ui appearing in this expression is an operator, which in the so-called einbein approximation 

dM 
du)i 

rsh m 



is defined by an extremum condition — 0. A simple expression for the spin-averaged mass M(nl) follows from the 



, ,. Wl LOo TO? TOo „ , . 

M{nl) = -l + -l + -^ + _L + E nl (u Ted ). (23) 

Z Z ZU1\ ZW2 

Here, the excitation energy E n i(uj Te d) depends on the reduced mass w re d = ^i+u 2 ■ The formula ([2"B")l does not contain 
any additive constant; for a light quark (e.g. in the D-meson) a negative (not small) nonperturbative self-energy term 
appears, proportional to (w u ) _1 ; it has to be added to their masses [28]. In the case of charmonium this term is small; 
the variables LOi(nl), the excitation energy .Enj (w re dK and the w.f. are calculated from the Hamiltonian (|21[) and two 
extremum conditions d M(nS)/du!i = (i = 1,2), [9ll2^|: 

H fnl(r) = M(nl) tp n i(r), 
2/ n 2 2ujf(nl)dE(nl,uj Ted ) 

^ {nl) = m * d^T) ' il = 1 > 2) - (24) 
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The potential Vb (t) in (|22l) is derived in the framework of the Field Correlator Method [Iol - [l2| and is the sum of a 
pure scalar confining term and a gluon-exchange part, 

V B {r)=ar- , (25 

A r 

where the vector coupling ck_b(Y) is taken in the two- loop approximation and possesses two important features: the 
asymptotic freedom behavior at small distances, defined by the QCD constant As(n/) [which is considered to be 
known, because is directly expressed via the QCD constant A-pg(ny) in the MS renormalization scheme]; it 
freezes at large distances. Details about the effective fine-structure constant can be found in Ref. [12[. 

III. RESONANCES IN THE DECAY SECTOR 

As it was discussed in @, the situation of two coupled sectors I, II, QQ and (Qq){Qq), can be treated in two ways: 

1) as a coupled system of matrix Green's functions, 

G ab , a, b = 1,11 

2) as a reduction of two-sector problem to the one-sector problem with energy-dependent "potential" V121 or V2i2 3 - 
We shall continue our one-channel treatment from the point of view of Sector II. In the same way as it was done 
before, one can define the "potential" V212 = Ki 2 n 3 ,n 2 n 3 

V mn3 , nWa (p, p', ^)=E J " + " 2 "i P ^ (P ° - (26) 

n 

Defining also J„„ 2 „ 3 (r) = J ^Jpr e 4 P r J„„ 2 „ 3 (p), one can write 

Vn 2 n 3 ,n' 2 n' 3 \r, F , - ^ E - E- 

n n 

and as a result one obtains a system of equations in sector II 

(r, r', E) i> n '-n'(* f ) d r' = (r) (28) 



where V22C 1 *) is a direct interaction between two color singlet mesons, which we neglect in the first approximation, 
Hq is the same as in (I2ip but for m,\ y m% equal to masses of mesons with quantum numbers n2,ri3,Ho = -ffo( n 2^3)- 
Neglecting V22, one can easily rewrite (|28|) for the separable interaction (|27|) 



= -£ / dVdVGW ,(r.O J ™™W J p<W Mn ,{r»), (29) 

n •* n 

where 

p j3i, „ik(r-r') 

(30) 



Introducing ip n = J J„„ 2 „ 3 (r)^„ 2 „ 3 (r) d 3 r, and integrating both sides of ([2Uf with J m „ 2 „ 3 (r) dr, one has from (|29l) 

!mn(E) <£ 

E — E n 



^ = 1^ (31) 



3 Note, that a parallel treatment of the open channel problem in nuclear reactions is developed by H.Feshbach with the help of the 
projection operators in his unified theory of nuclear reactions [3dl |. 



7 



with the same w mn (E) as in ([8]), and the equation for eigenvalues is again (JTTJ) . 

Since the CC interaction (|27p is separable, one can study the structure of the spectrum of our CC problem in 
more detail; in particular, whether there can appear poles (CC resonances in terminology of |5||) due to strong CC 
interaction, which are additional to the one-sector spectrum of poles E n , the latter being simply shifted by CC. As 
was argued in [5[, we define the integral 



\J n 



,(P)| : 



W^rT (E-E n )(H^\p)-E) 



(32) 



According to a bound state in a single channel 7i2?i3 due to CC with the sector I can exist, if in the region, where 
is real (below threshold E t h{n2,n^j), it becomes larger than one 



I n2 „ 3 (E)>l, E < E th (n 2 ,n 3 ) 



In the momentum space one has 

H (p) ipn 2 n 3 (p) 



,n'n'AP, P', E) 1p„>„>(p'] 



d 3 p ' 

(27T) 



Elpn 2 n 3 (p), 



(33) 



(34) 



which yields the same equation as in (I3T1) . As before in Eq. (Till) , one obtains from ([34]) the equation det \E—E n — w) = 
0, which defines all poles in the cut E plane below thresholds and on the second, and higher Riemann sheets. 



IV. THEORY OF COUPLED-CHANNEL RESONANCES BASED ON THE WEINBERG EIGENVALUE 

METHOD 

The coupled-channel problem can be quantified using the eigenvalue analysis introduced by Weinberg 20] . Although 
this formalism has been developed long ago, it is still not widely known. That is why in this section we present a 
short summary of corresponding formulae leaving details to the Appendix D. 

The Schrodinger equation for two-body like Hamiltonian H = Hq + V can be written in the standard (time- 
independent) way 

(H -E)V E (r) =-V* B (r), (35) 

where E is a spectral variable, Vl/£;(r) is an energy eigenstate and V is an operator, which in the nonlocal case acts in 
([33]) as ^(^^(r)) = J V(r, r')4 , £;(r') dr' . In the Weinberg method, instead, w.f. are the eigensolutions of 

(H -E)^ v {r,E) = ^-* v {r,E) (36) 
r)v{E) 

where E is the continuous parameter entering w.f. and the index v labels the discrete eigenvalues and eigenvectors. 
The Weinberg eigenvalue r\ v {E) is the potential scale and thus the spectrum consists of all the potential rescalings 
that give solution to that equation, for given energy E. 

Let us now turn to the question of the rigorous definition of the resonance wave function and start with the one-state 
situation, when only one state is considered in sector II, with fixed n2,"3- The induced interaction V2i2(r, r', E) has 
the form (|2"T1) and direct interaction V22 is neglected for simplicity. One can exploit the WEM [2(|, and divide the 
potential V2i2(r, r', E) in (|32|) by an energy-dependent factor rj u (E) considering instead of Eq. (|28|) another one (with 
V22 = 0) 

H *„(r, E) + [ V ^ V ^ E ) y v{r > E) d 3 Y > = E $j/(r E ) , (37) 
J VA E ) 

which defines for each E the eigenvalue r) u (E) and eigenfunction ^ „{r, E) 1 with the boundary conditions 

<M0) = const, ^(r -> 00) = C - — , k = ^2M(E - E th ). (38) 

For E < E^ one has instead ^ u (r — > 00) ~ c exp(— nr)/r. In the WEM, resonance structures as well as bound 
sates can be obtained in terms of Weinberg eigenvalues r] v (E). Note, that a solution of integro-differential equation 
(|37p in the coordinate space can satisfy these boundary conditions for each energy E only for some discrete value of 
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^(E), v = 1, 2, ... . Compare e.g. with the case of bound states (energy below threshold), where boundary conditions 
at origin and infinity can be matched only for discrete energy E^ in standard formalism (rj = 1) or at some rji(E) for 
any E in WEM, with the relation rji(Ei) = 1. 
The normalization of wave functions is 

J drdr'* u {T,E)V 2l2 {T > r',E)^ u ,{r,E) = -5 vu ,7 lv (E), (39) 

Note, that ^212(25) is real analytic (holomorphic) for all E except for pole positions, and the off-shell i-matrix looks 
as (see [20j and Appendix D for a derivation) 

t(p, p', 25) = - £ = 71 AE ^, <v(p, E) a„(p', E), (40) 

with 

a„(p, E) = (H (p) - E) *„(p, E), (41) 



a„(p, E)a v >(p,E) d 3 p _ 

H (p)-E (27r)3-^' ! ^ 



and the Green's function (JTDJ) has the form 



^W 1 - 2 > s ) - X — Y^ae) — • (43) 

The sum over v is fast converging as one can see from the example of square well and other potentials fast decreasing 
at oo 21]. Therefore in what follows in our calculations we shall consider only one term in the sum over v, which is 
relevant for a given threshold. 

The Breit-Wigner resonances in sector II are obtained from the condition that for some v = vq, r\ VQ {Eq — ^) = 1, 
and 

Vuo (E) = l + rl va (Eq - y ) (e - E + + ... (44) 

Note, that the corresponding (r, E) serves as the normalized resonance wave function and can be used e.g. to 
calculate average values of some operator or perturbative shift of resonance position. The QQ Green's function with 
account of channel coupling can be written near the pole E = E R (resonance) as 

(z) _ *,(!,£) *+(2,£) _ tt„ 

qq( ' ' > ~ e ' *»- / dvAER) ^ 

V dE 

As it is seen from (13"T|) . the introduction of WEM eigenvalue in equations reduces to the replacement w mn (E) — > 
w ™?Ja , hence the resulting equation for the calculation of rj(E) is 

det (l - ™ ( f ) 1 A = 0, (E) mn = E n 6 mn . (46) 
V V{E)E-Ej 

Eq. (l4^)l is of the n-th power in r], when n levels (QQ) are taken into account, and this yields n roots r]k{E), k = 1, ...n. 
Total number of poles is given by solutions r)k(Ei) = 1, I = 1, ...N p , where -/V p depends on behavior of r)k(E). 

We started formally with the one channel in sector II, i.e. with fixed ri2,ri3, hence in w(E) in Eq.([4"6"| the sum 
over ri2,n 3 (cf. Eq.®) reduces to one term. However, for several states n 2 ,n 3 one has equation (f28f with interaction 
kernel V212 as a matrix in indices 7^2,^3, n 2 ,n' 3 , and if Tj(E) in (|3T[) does not depend on 712,713, then as a result one 
has for rj(E) the same Eq. (|4"6"|) . but now with w(E), which corresponds fully to ©, i.e. contains the sum over 7i 2 ,7i 3 . 

Let us discuss how the basic equation (|4"o| changes for many channels ni , 773 . We start with Eq. (j2TJ)) , which is 
equivalent to (|3"T)) . when one introduces in (|29p in the denominator on the r.h.s. the factor r\ v (E). Multiplying both 
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(49) 



sides of this modified Eq. (|29[) with J m n 2 n 3 (r) and integrating and summing over 722,713 one obtains equation similar 
to ([51) 

' (47) 
where (p m (E) = Jmn 2 n 3 (r) ipn 2 n 3 (r) d 3 r, and w mn {E) is the same, as in ([5]), i.e. again with the sum over 212,213. 

The resulting equation to determine ij(E) is again (|46|). and all equations (l3"7T) - ([3"9")) have the same form, if one takes 
into account, that ^> v is a column of ipn 2 n 3 components and V212 is a matrix in indices 72271,3, 72 2 ?2 3 . 
Finally, the separate components ip n2 n 3 are found through ip n (E) via (cf. (I29p ) 

^n 8 (r) = - E vA Ef(E-E n ) I ( r ' r ') (*') (48) 

and partial widths of the resonance are found in lowest approximation as (for one channel 72 in sector I) 

T nn2 n 3 (E R ) = 2Im n2n3 (w nn (E R )) = 2ttJ |J„„ 2 n 3 (p)l 2 - En 2 n 3 ( P )). 

To understand the possible origin and position of resonances in our CC problems, one can consider several typical 
cases, depending on relative positions of bare resonances E n and thresholds E t h (722223). Consider first one state in 
sector I, one state in sector II, then w nn {E) < for E < E t h, and Re (w nn {E)) changes sign at E = E* . The resulting 

qualitative picture of r](E) = w £"^ is shown in Figure [2] for three cases: E n < E t h (Fig. EJa)); E th < E n < E* (Fig. 
[2jb,d)); E n > E* (Figl2](c)). In Fig[2]Ja,b) one can see one critical energy for which rj(E) = 1. This point corresponds 
to the shifted energy level E n . 

Note the possibility of a pair of additional roots of equation 77 (E) = 1, when E n > E* , E t h < E n < E* and 

w nn {E th ) > \E n - E th \. (50) 

The condition ([50]) defines the strength of CC interaction in the situation depicted in Figl2Jc,d) which is necessary 
one additional pole near the threshold energy (see Appendix E for details) . As we shall show below in section IVI1 the 
situation of Figj2](d) is most likely realized in the 3 Pi state of charmonium, where the threshold peak corresponds 
to the A (3872) state. In this case actually two close- by thresholds are present (DqDq and D + D*_), and as will be 
seen, the experimentally observed situation with one peak at lowest threshold and wide structure near E ~ 3.940 
GeV indeed occurs. 

Consider now the case of two levels in sector I, E\ and E2; and one (or more) state in sector II. The equation for 
77 (E) has the form 

rj 2 {E) - rj(E) (w u + W22) ~ *12 w 2 i = 0, (51) 



with wik = 'wzf i and the result 



V± i E ) = ^ (*11 + W22) ± ^ V (Wll - 7«22) 2 + 477J12 W 2 1 ■ (52) 



where for notational convenience we have suppressed the energy dependence of the w. 

E-Ex t *" u 'W^J — E-E 2 



Near E = E%, rj±(E) can be identified with the one-channel eigenvalues rj\{E) = and 7/2 (E) — W22 ^ 



namely for E < E\ and E — > E\ one has 

it? , t? \ ( J7 l \ _i w 12 {E)w2i{E) 
w n (E) (E - E 2 ) 

^ $20% (54, 

and for E ->• E 2 one should change in ([53]), ([M]) 1^2. 

The situation with trajectories rj± (E) is in general rather complicated, and we describe below only one case when 
Ei < E 2 < E* k ,i, k = 1, 2 where Re (wik(E* k )) — 0, and in case of strong mixing of channels 1 and 2 the point Eq, 
where Re (77+ (-Eo)) = Re (r]-(Eo)) lies between E\ and E 2 - (The position of E t h is irrelevant for the situation where 
all imaginary parts are neglected). However for weak mixing of channels 1,2 roots rj + ,r)_ never coincide. One can 
see from (|52[) that in the weak mixing case only two poles remain, corresponding to shifted levels E\ , E2 and no new 
resonances appear at least for E < E* k . 
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(c) case E n > E* (d) case E th < E n < E* 



FIG. 2: Qualitative pictures of Weinberg eigenvalues Re(r;(_E)) as a function of energy E, where E n is the eigenvalue of the 
single-channel relativistic string Hamiltonian Ho (bare resonance), E t h denotes threshold and E* is the point, where Ke(w n n(E)) 
changes sign. 



V. MIXING OF STATES IN THE WEINBERG FORMALISM 



The WEM solves the important problem of constructing the full set of orthogonal states in the coupled channel 
problem, and thus the problem of mixing of states. This is nontrivial in the situation under investigation, since the 
interaction in the sector I induced by the coupling to the sector II, Vi 2 i(r,r'), Eq. ([TBI) , is energy dependent and 
hence violates the orthogonality of eigenstates. In addition, for energies above threshold, this interaction is complex 
and makes the corresponding states the resonances, which cannot be normalized and orthogonalized to each other in 
the ordinary way. Happily, WEM allows to define all states and their mixing in the mathematically rigorous way, as 
we shall now show. 

We start with the formulation in sector I and write starting from (P2U1) the WEM equation 

H Q vMr, E) + [ Vl2l{r ) Y ': E) <Mr', E) d 3 r' = E ^(r, E) , (55) 
while the unperturbed states ^ n (r) satisfy 

# *„(r) = E n 9 n (r). (56) 
Note, that ^ v {v,E) depend on energy E, while vE f n (r) do not. Similarly to ([39)) . the orthogonality condition is 

f drdr'y„(r,E)V 1 2i(r,r',E)y u ,(r',E) = -8„„, Vu {E). (57) 

Consider now the expansion of a WEM state in the set of \l/ ra states, 

#„(r, £) = £<(£)*„(r). (58) 

n 
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Taking into account, that 

J drdr'y n (r)Vi2i(r,r',E)y m (r')=w nm (E), (59) 
and multiplaying both sides of (|55|) with \lv (r, E) and integrating over dr, one obtains 

c n ( E )c n (E) {E n -E) +2_^ —t^t =0. (60) 

n m,n ^ ' 

Thus one obtains the equation for eigenvalues tj u {E) 

dct(£-£ + -y= (6i) 

which coincides with (j46|) . obtained in sector II. Now we are specifically interested in the coefficients {c^}, {c^ } for 
two different eigenvalues •q u (E) 1 r) v i(E). 
The first condition follows from ([57)) . ([60]) 



Y J <l{E)C{E){E n -E)=8 vv ,. (62) 



It is convenient to introduce reduced coefficients: 



c - (i?)= 7^' %n(£)= vFin - (63) 



Then the solution for two eigenvalues in (fBTj) is 



/■T-.S -((w 2 2 + wn) ± ^(w 2 2 + wn) 2 - 4det w) 
«M#) = — 2 *- ' 

here and further for notational convenience we will suppress the energy dependence of the w nm . The normalization 
condition has the form 

J2 K(E) ci{E) = <W, J2 <W ™nm(E) c" m {E) = -«W (65) 

n n.m 

Let us take one concrete example of two states in the subthreshold region (e.g. (2 3 Si) and (l 3 Di) states of charmo- 
nium, however at this stage they are not specified). 

Keeping only two states n= 1,2 e.g. for (2 3 Si) and (l 3 Di), one can write for c^(E), v = a, j3 

c*(E) = cos(p{E), c*{E) = smip{E); c((E) = ahxip(E), c 2 {E) =-coacp(E). (66) 

Note, that the appearance of 0(2) coefficients is not accidental since w nm (E) is symmetric in n,m. 

We are thus e.g. looking for the shifted and mixed (2 3 Si) state, denoted by a, and the same for (l 3 Di) state, 
denoted by (3. 

,.« n _ «»y(*Q ^ , ^(E) , T ,_ 

V E\ — E y — ^ 

s, . sin<z>(i?) cos <p(E) 

To find cosip(E), one can use the second equation in (|65p. which yields 

sin 2 </?(-E) wn — cosip(E) (w 12 + w 2 i) + cos 2 y{E) W22 — —1113(E) 



cos 2 



(p(E)wn - cos (p(E) (w 12 +W21) + sin 2 tp(E) w 22 = -T) a (E). (68) 
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This gives the condition W11+W22 = — { r lct(E)+r]p(Ey\, which is identically satisfied, and the final result for cos 2 <p(E) 

tin - w 2 2 + D 



cos 2 ip(E) = — , D = \J (wn - W22) 2 + 4wi2W 2 i- (69) 

Note, that the sign of D is connected with the corresponding choice of the root in for rj a (E) (lower in energy 
state) we have chosen the sign +. 

It is clear, that cost/? depends on E and therefore to define finally the mixing coefficient, one should fix the energy. 
E.g. for the state (3, the eigenvalue rjp{E) crosses the line r](E) = 1 at the resonance position E = EE, complex in 

general, and the mixing coefficient of interest from (|66|) is cP x — s ™ v ( e p ( while the mixing coefficient of the state a 
is to be taken at E = E%, c? = si "' p( ' E ° ) . 

Hence, for small shifts EE = E2,E* ~ E\, and energy independent ip, one recovers the symmetry condition 

|c?|«|c£|. (70) 

Finally, one should connect normalizations of ^ n and y& a 'P . This can be done, if one considers the limiting case of 
one channel v, where according to (|5^|) . (j57| . one has 

(C(E) w nm (E) c» m {E)) = - Vu (E) (71) 

and for E = E R (at the resonance position), T] V (E R ) = 1, and for one level n from (f61j) one has w nm {E„) = E R — E n = 
—AE n . Hence in the one-channel - one- level limit we have 

«) 2 A£„ = 1, < = ^=L=. (72) 
Therefore if only one level n is kept, then the normalized WEM states can be defined as 

* a (£?«) = (E* ) y/AEZ, J (* Q (E* ■)) 2 d 3 r = 1 (73) 
and finally the standard normalized mixing coefficients are 



E 1 -E* p ' 2 \/E 2 - E* 

where AEp — E2 — E^, AE a = E% — E R . One can see, that in general coefficients are less than unity due to ratios 
of square roots. We finally write for sin ip 

(wn - W22) 2 + 2 (w 12 ) 2 + 2 (w 2 i) 2 - (ti>n ~ W22) D 
2D' 



sin 2 p(E) = \ ^ W2ZI Z^l^l ^ W22J ^ \ (75) 



Another (and physically more motivated) normalization for ^"(i?*) follows from (|45|) . which can be written as 

G ( nUl,2;E) 



J QQ^^^> r E R _ E \ indE^l 



dE 



Estimating dWn ^ ^ — j§^r^j, £ < 1, one obtains dVa } E a ^ — E ]_ E , + \E R -E th \ > anc ^ tne defacto wave functions is 
^ a {l,E R )/J dVa ^ R \ which is close to for £ < 1. 



VI. RESULTS AND DISCUSSION 



The formalism given in this paper is based on the explicit knowledge of wave functions in both sectors I and II 
and yields the CC interaction operator w(E) expressed via the overlap integrals, see Eq.©. The resulting effective 
interaction in each sector is energy dependent due to w(E), and violates usual orthonormality properties for wave 
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TABLE I: Charmonium spectrum in the single-channel approach derived from RSH (|2ip [T3 |. The experimental numbers are 
taken from PDG [3l[. All masses are in GeV. 



State (Thresholds) Theory Experiment 



IS 


3.068 


3.068 


IP 


3.488 


3.525 


25' 


3.678 


3.674 


(DD) 




3.729 


ID 


3.787 


3.771 


(D*D) 




3.872 


2P 


3.954 


3.930 


(D*D*) 




4.014 


3S 


4.116 


4.040 



DD 




p, GeV 

FIG. 4: The Weinberg eigenvalue Re (77(E)) for 3 3 5i 
FIG. 3: The squared overlap integral § J2 I -Annans (p)\ 2 state, where £ n =4.116 GeV is the eigenvalue of the single- 

v* channel relativistic string Hamiltonian Ho (bare reso- 

for 3 3 5"i state. nance) and E t h(D* D*)=4.014 GeV denotes the closest 

threshold. 

functions. Moreover, new states appear for energies above thresholds, and one needs a rigorous formalism to treat 
the complete set of eigenfunctions for such operators. The WEM is indispensable for this purpose. In Eq. (l46|) explicit 
conditions are written down for Weinberg eigenvalues r)(E). It is important, that rj(E) has simple analytic properties 
in the S-plane. Therefore physical quantities expressed via tj(E), like scattering amplitude (IE.1|) or production cross- 
section (|E.2|) have a definite analytic expression near the pole(s), different from the Breit-Wigner form in general. 
This property is more important in case of the complicated arrangement of thresholds and poles, as it is in the case 
of X (3872), see below. 

Another practical advantage of WEM is the complete set of states for each energy E, allowing to define unambigu- 
ously symmetric mixing coefficients, as it was explained in section fVl 

Before a detailed discussion of results, one should stress two main features of the closed channel pole E n behavior 
under the influence of CC: 1) CC is attractive for all states below CC threshold; 2) CC is attractive in some region 
Eth < E n < E* above threshold and repulsive for E > E* (in the limit of small width) . Both statements follow from 
the condition Re (w(E)j < or Re (w(E)) > in ©. As will be seen in the simplest case of one channel with lowest 
threshold, the 2 3 P\ pole E n occurs in the attractive zone of the DD* channel and hence moves down with increasing 
coupling. 

Below we give several examples of WEM application to different problems in CC dynamics. We shall consider 

i) How CC interaction changes n 3 Si states as compared to one-channel calculations. We will calculate energy shifts 
and widths for 3 3 <Si state and also mixing between 3 3 Si and 2 3 Si states. 

ii) We calculate eigenvalues and amplitudes in the 1 ++ state in connection with the bare 2 3 P\ level and resulting 
AT (3872) resonance. 



14 



TABLE II: The decay width T (p) of Z^Sx char monium state. The resonance momentum p is taken from PDG [3lj . 



Channel p, GeV F(p), MeV 





DD 0.777 0.31 
D*D 0.576 25.5 
D*D* 0.227 17.8 






TABLE III: Ratios of branching fractions for 3 3 Si state. 




Ratio 


Experiment [32J This paper 




23(^(4040) 
B(ip(4040) 


DD)/B [ip (4040) ^ D*D) 0.24 ±0.17 0.012 
D*D*)/B(ip(4040) -> D*D) 0.18 ± 0.18 0.70 


0.003 
1.0 



To illustrate this formalism we will consider situation with one level in sector I and one (or many) level(s) in sector 
II. In Table U we present charmonium mass spectrum in the single- channel approach (SCA) derived from RSH (|2ip 
(see for example |12j |) in comparison with experimental data and showing the thresholds. 

A. 3 Si levels 

As a first numerical example we consider the mass shifts and widths of the n 3 Si (n — 1,2,3) states. For these 
levels the corresponding y\g factors are y r ilt( 3 Si -> DD) = ft§j}( 3 Si -> D* D) = i e ijm q m ; y{^S 1 -> D* D*) = 

{pij Ik ~ fijk 1i + $ik Qj) ( see Appendix B,C) and the transition matrix element is rewritten in the following form 
(see appendix C and equation (|C.2I0 

J ni „ 2n3 (p) = -2L| ^y 1 g(p,q)^ ) ( C p + q)^ ) (q)4 , ; 3) (q), (76) 

where 7 rj 1.4 is channel coupling parameter which is proportional M u (see Appendix C) and c = - ^ ~ 0.73, 
where the averaged kinetic energies of heavy and light quarks in D meson u> q ~ 0.55 GeV, wq ~ 1.5 GeV are taken 
from In Eq.fTBJ tf^, ... = /Vilr, ... are series of oscillator functions, which are fitted to realistic w.f. (see 
Appendix A). We obtain the latter from the solution of RSH (HD (Hj]. 

The widths and mass shifts are obtained from \ J ni n 2 n 3 (p)| 2 averaging over initial (i) and summing over final (k,j) 
polarizations. Note that the final formulas for the width in channels DD, DD* and D*D* differ by spin factors, 
which yield the ratio 1:4:7. From Eq. (ll3[ I4U)) one can write the width taking into account relativistic corrections 

F ni n 2 n 3 (P) = ~ l^m 2 n 3 (p)| 2 f — F=f==r + — F=f==r J > (77) 
n \ VP + M n 2 y/P + M nJ 

where M„ 2 , M„ 3 are the masses of the corresponding D mesons. It is important that the value of the decay width 

TABLE IV: Mass shifts (in MeV) of the n 3 Si states with n = 1, 2, 3. 



State DD 


DD* 


d*d* 


Total 


l 3 Si -5 


-19 


-30 


-54 


2 3 Si -15 


-41 


-56 


-112 


3 3 Si -6 


-10 


-45 


-61 



strongly depends on transition matrix element. This is illustrated by the behavior of \ J ni n 2 n 3 (p)| 2 for state. 
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(a) One threshold, E th (D D* )=3.872 GeV. ( b ) Two thr esholds, E th {D Q D* Q ; D+D*_ )=3.872; 3.879 

Oe V . 

FIG. 5: The Weinberg eigenvalue Re (r/(E)) for 2 3 Pi state with different values of channel coupling parameter ((1) - 7 = 0.6, 
(2) - 7 = 1.0, (3) - 7 = 1.1, (4) - 7 = 1.2, (5) - 7 = 1.3), i? n =3.954 GeV is the eigenvalue of the single-channel relativistic string 
Hamiltonian Ho (bare resonance) and E t h(DoD*,; D + D*_ )=3.872; 3.879 GeV denote thresholds. 



As it can be seen from Figure [31 \Jmn 2 n 3 (p)| 2 is oscillating and has two zeros, corresponding to the wave function 
nodes. In the small width approximation (TT31 the width and shift of E n level will vanish when p (E n ) approaches zero 
on Figure [3] It is not a physical situation, and in the next approximation one should solve Eq.([3j) in complex plane 
and take into account possible mixing between states due to open channels. For instance it can be 3S-2S, or 3S-2D 
mixing. Due to the mixing, the w.f. of the "pure" states changes and minima in Figure [3] can be filled in by admixed 
states. In Tables [TTl IIIII are given the small width values for 3S state of charmonium in the D D channel, illustrating 
the zeros discussed above. 

In WEM the shifted level positions are defined from Eq. (14T>)) and for 3 3 Si one obtains the picture shown in Figure 
SJ The level shifts calculated from Eq. (fT3")) are given in Table ITVl One can note relatively small shifts (AE < 100 
MeV) as compared to [f| 0], where 3 Pq and SHO model was used, whereas in our case more complicated realistic wave 
functions were exploited. 

In addition we have considered mixing between 3 3 Si and 2 3 Si levels via D*D* threshold, which turned out to be 
small, with the mixing angle (defined as in ([7S]) ) ip = 5°. 



B. 2 i P 1 -level 



A separate discussion is needed for the s-wave decay to charmed mesons. We take as an explicit example the decay 
2 3 Pi — > DD*. Note, that due to positive C parity the s-wave strength is mostly concentrated in the DD* channel. In 
this case, the situation of Figure [UJc,d) is realized when Re (r)(E)) can cross the unity line at several energy values, 
thus producing several resonances. In our calculations we show Re (rj(E)^J in Figure [5] which correspond to different 
values of channel coupling parameter in the region ±30% around the standard value 7 = 1.4 (M u = 0.8 GeV). As it 
can be seen, Re (r](E)) intercepts the line Re (j](E)) = 1 three times. However we have to take into account imaginary 

parts above the thresholds. The simplest way is to calculate factor rjj^gjp which appears in the squared t-matrix 
(|40|) . The result is the two-resonance structure, one of which is near threshold M ~ 3.872 GeV and another one near 
M ~ 3.940 GeV, the latter becomes increasingly broad with increasing coupling 7 to open channel. In the recent work 

33] a similar form of the first peak was suggested. We note, that the factor rjz^rgyp is relevant for the t-matrix of 

DD* scattering, while new charmonium resonances were observed in production cross sections like e + e~ — > DD* or 
B —> K(DD*). Therefore we define the production yield ^(E 1 )! 2 , given in (IE.2[) and show in Figure |5]the quantity 

\E-^E lW -'w^ E (E)[ j! ~ \i-ri(E)fi(E-E ) ' ^ n our approximation (DqDq and D + D*_ thresholds coincide and there is no 
connection to uJ/^/) and J/'ipirir channels) one can see the double peak structure for 7 = 1.0; the first peak at 3.872 
GeV is accompanied by a wide peak around 3.940 GeV. However, with increasing 7, when 7 = 1,2, the peak in 
Fig. [S]at 3940 becomes flat, while the lower peak at 3.872 GeV is narrow and high. This picture corresponds to the 
experimental situation. 

The case when both thresholds D Dq and D + D*_ are taken into account, is illustrated by Figure[S]for Re (^(P)) and 
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FIG. 6: Production crossection ~ i^,™™^"^^^ for 1 ++ state with different values of channel coupling parameter ((1) - 
7 = 0.6, (2) - 7 = 1.0, (3) - 7 = 1.1, (4) - 7 = 1.2, (5) - 7 = 1.3). For small values of channel coupling parameter 7 (line (1)) 
one can see a good Breit-Wigner shape, which corresponds to the shifted 2 3 Pi state, while for larger 7 (lines (2), (3), (4)) there 
is a broadening of higher resonance together with steep rise near the threshold Eth(DD*) = 3.872 GeV. 



Figure [7] for the production cross section. As can be seen, the curves for production cross section depend strongly on 
channel coupling parameter 7. For 7 = 1.0 (line (2)) which is 30% smaller than the nominal value 7 = 1.4 (M u = 0.8 
GeV) one can see a peak at the higher threshold D + D*_ and a wider peak at 3.940 GeV, however for 7 = 1.2 (line 
(4)), the 3.940 GeV peak flattens and simultaneously^), the peak appears at the lower threshold DqDq, while only 
a week cusp is seen at the higher threshold D + D*_. Surprisingly, the isotopically equivalent thresholds (which we 
take into account with equal weight) due to different position in energy plane, provide finally the asymmetric picture 
observed in experiment [221 ] . 



VII. SUMMARY 



We have formulated equations for Green's functions of strongly coupled sectors, where new resonances can appear 
due to CC interaction. We found that the best formalism for the CC induced energy-dependent interaction is the 
Weinberg eigenvalue method. Conditions for the poles and their positions were systematically studied in case of 
P-wave and 5*-wave channel coupling. In the first case one finds only displacement of poles, while in the second new 
resonances appear, and in the 3 Pi case two peaks at 3.872 and 3.940 GeV were found with the height depending 
on the coupling constant M w . Moreover, we have shown in FiglZl that at one value of M w the lower peak is at the 
DqDq threshold (but not at the D + D*_ threshold) and at the same time the upper peak at 3.940 GeV flattens. This 
situation corresponds to the experimental data [22| and supports our dynamical CC mechanism. 

Mixing of n 3 Si states was formulated in WEM and found to be small, while shifts of 3 3 Si are of the order 50-80 
MeV, which signals necessity of mass renormalization. 

The method developed in the present paper, provides a rigorous definition of resonance wave functions and mixings 
in the case of strongly coupled channels. 
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FIG. 7: Production cross-section ~ [Birp^zSf^fgTp ^ or ^ ++ s ^ a ^ e with different values of channel coupling parameter ((1) - 
7 = 0.6, (2) - 7 = 1.0, (3) - 7 = 1.1, (4) - 7 = 1.2, (5) - 7 = 1.3) in case, when both thresholds DqDq and D + D*_ axe taken 
into account separately. For small values of channel coupling parameter 7 (line (1)) one can see a good Breit-Wigner shape, 
which corresponds to the shifted 2 3 Pi state, while for larger 7 (line (2)) there is a broadening together with a cusp first near 
the closest threshold EthiD+Dl) = 3.879 GeV and then for 7 = 1.2 (line (4)) a sharp peak appears at the Eth(DoDo) — 3.872 
GeV. 
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Appendix A 
WAVE FUNCTIONS 



In Eq.ra R ( "X\ Rn? ] and R { ™ 3) are series of oscillator wave functions, which arc fitted to realistic wave functions. 
We obtain them from the solution of the Relativistic String Hamiltonian (EH]) , described in [9|, |12|, |25( . 
In position space the basic SHO radial wave function is given by 



Rtf°(fi, r) = P'^ ^ltllfi W? L n-l\P 2 r 2 ) (A.l) 

CO 

J (R s n f°(l3, r)f 



r dr = 1 



where /3 is the SHO wave function parameter, and L l ^( 2 (j3 2 r 2 ) is an associated Laguerre polynomial. The realistic 
radial wave function can be represented as an expansion in the full set of oscillator radial functions: 
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TABLE V: Effective values ft (in GeV) and coefficients Ck of the series of oscillator radial wave functions i?f i ff °(/3, r) which are 
fitted to realistic radial wave functions R„i(r) of charmonium and D meson. 

State /3 Coefficients Ck 

Charmonium 

IS 0.70 ci = 0.97796 c 2 = 0.169169 c 3 = 0.117682 c 4 = 0.019694 c 5 = 0.025113 

2S 0.53 ci = -0.11889 c 2 = -0.972774 c 3 = -0.134041 c 4 = -0.142303 c 5 = 0.000142 

35 0.46 ci = -0.09354 c 2 = 0.149573 c 3 = 0.958816 c 4 = 0.112102 c 5 = 0.183886 

2P 0.48 ci = -0.06271 c 2 = 0.981834 c 3 = -0.123392 c 4 = 0.127111 c 5 = 0.000495 

D meson 

IS 0.48 c=l 




= £c fcJ R£f°(/?,r). (A.2) 
fe=i 

Effective values of oscillator parameters j3 and coefficients Ck are obtained minimizing \ 2 an( i listed in the Table [V] 
The quality of approximations can be seen from the Figure HI In the momentum space the SHO radial wave function 
is given by: 

R nl IP. P) ~ ^575 

OO 





/ 2(n-l)! 
T(n + I + 1/2) 



V/2/3 2 ri+1/2 



Appendix B 
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THE VERTEX OPERATORS AND SPINOR EXPRESSIONS IN THE (2 x 2) FORM 

Our purpose here is to go from Eq.©, where y\fi is the trace of (4 x 4) form to the (2 x 2) or spinor form, defining 
in this way y^i- 

_ <-> 
We consider operators of the form (■0r i '0), with Tj consisting of Dirac matrices 7^ and derivatives di- To proceed 

to the (2 x 2) form, one exploits the limit M — > 00 of the heavy quark mass, so that for the light quark in the 

heavy-light meson the Dirac equation can be used, and one can use symbolically Dirac one-body equation, (ap + 




(3(m + U)) ip = (e — V) Vj so that for ifi = I J , one has w = m+u 1 _ v+e (<rp) v. One can also use connection 



ip = C- 1 tp c = ip c {C' l ) T , where C = {C' l ) T = 7274, 7* = -i/3oti, so that ipTiip = (v c , w c ) 7274 I I- Note ' 
TABLE VI: Bilinear operators tpTiip and their (2 x 2) forms (Notations see in the text). 



jpc 2S+i Lj ^ (2 x 2) form. 

^ So —ijs v c v — w c w 

1" 3 Sl 7i — (v C (TiV + W C (TiW) 

1+- 1 P 1 -ijsdi v c *d l v~w c 'd i w 

0++ 3 P 1 i{v c w-w c v) 

1 ++ 3 Pi 7i75 -(v c a l w + w c a x v) 

2 ++ 3 P> 7i dfe +7fe 9i -^S ik B -(v c p ik v + w c p ik w) 



2"+ 1 d 2 (didk-kSik(d) 2 h. 



75 



2 " i D 2 (rti dk +lk di -\5ikd)^ -(v c p lk w + w c p ik v) 
1 3 Pi 'jiOJik ~(v c aiLJikV + w c aiOJikw) 



that spin indices of charge-conjugated spinors are connected to ordinary spinors by matrix 02: v c 02 — — (02 v c ) T = v c , 
and for w c one has: 

w c <r 2 = ( 77— r r crpv c ) o 2 = -v c cr^p 1 = -w c (B.l) 

\m + U-V + s J m+U-V+e y ' 

where the notation "Jf implies, that operator acts on the left. We are considering 7 lowest states and display in 
the Table IVT1 the operator T.j, the corresponding quantum numbers J PC , spectroscopic notation 2S+1 Lj and the 
equivalent (2 x 2) form for the same vertex I\ in the last column. Wc are using in the Table IVTl the following notations 

pik = &i dk +&k di di Sik 4 , Lo%k = didk — ^ S lk (d) , d=di It- 
Note, that in the 2x2 form one has: 



m+U-V+e 



3 „ . _c Pfc C( 

Pi : -ie lki v — — — — 

m + U — V + s 

<-> 

D 2 : v — — — — uiikv 

m + U — V + e 



In the (2 x 2) form one can write wave function of charmonium and D-mesons (^qq , ipQ^ , ^q^) as 
Const ip n i(r) (y c T^ d and normalize it as 



I*, 



= 1 = | \^r)\ 2 r 2 drtr{T%yS}dV. (B.2) 
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TABLE VII: The V^J operator and nonrelativistic form yl%$ for DD, DD* and D* D* channels. 



jPC 2S+l Lj 


p( n l) 

red 




ylii 








DD 


-j= {DD* ± DD*) 


D*D* 



++ 3 -Po 72°" n TS^ 11 ) " 75 - 5jfc(qn)J 

1 ++ 3 -Pi ^ £ikio- k ni - i J§ Uinj - (qn)d«) 



2 ++ 3 P2 § (o^n; + oirn - § (<rn)<5iij - - j \ni{qkSij - qi5 jk + qjS ik ) 

+rii{q k 8ij — qiSjk + qj$lk) 
— n j(| gfe^ii) — rifc(| tfofo) 
+ |(qn)(5ii(5 ji fej 



Defining (p„i{r) = R n i(r)/y/An, the normalization condition for angular part takes the form J tr |r|,™^ r|," rf + j ^ = 1. 
Then Jnin 2 n 3 (p) can be written as in ([7]), but y^i can be found in spinor (2 x 2) form as 

S3( = to{rSrSHrH} (b.s) 

see Table IVlTl and all are normalized as written above. Hence e.g. 

^("2,3) / n \ _ ' r'" 2 - 3 Vnn - ^i. 

1 red — fct 1 red V 1 ^ 7 — 



Appendix C 
THE PAIR-CREATION VERTEX 



In the same way we consider here the (2 x 2) reduction of the pair-creation vertex, taking tp, tp for light quarks as 
solutions of Dirac equation and writing the effective string-breaking Lagrangian as 

C sb = [ ${u) M u <P(u) d 4 u = M u [ W f Pv d 4 u (C.l) 
J Jm + U-V + eo 

and we have denoted: U = or; V = — |^f-; P= p— P; Eq is the Dirac eigenvalue £o = M {Qq) —Mq, v c — v c oi is 
the spinor of antiquark, M u is the same as in Eq. ([3]). 

One can take in (|C.1[) the averaged value of the denominator (to + U — V + £o) m + (f ) — (V) + £o which 
effectively redefines our vertex constant M u . As a result the reduced form of the matrix element J{p) in Eq. ((HJ 
assumes the form 

J(p) = ^/ ^r^qJ^fcp + qJ^CqJ^W ( c - 2 ) 

where y\?£ is given in Table IVTA and 7 = m+ ^^y j+rg ■ One can find values of ([/), (V), £0 in Table lylTfl and 

persuade oneself, that 7 is rather stable for different a s and numerically 7 = 2 ®2G6v s=s 1.4. 

To check consistency of our approximation of putting average values into denominator, we have compared norm- 

lization conditions of bispinors {v + v) + (w + w) = 1 = v + (l + ^ e+ ^F^v)+my i ) u ' an< ^ found that the term with 
denominator contributes around 20%, and we expect the same accuracy in definition of 7. Actually, we are always 
varying 7 in the region ±30% around the nominal value 7 = 1.4. 

One can also check at this point how the (4 x 4) vertex M w j/123 goes over into the reduced form 7 y^i- E.g. for the 
1 state decaying into DD* one has in the heavy quark mass limit (see e.g. [13,14]). 2/123 = tg "J ifc " , and u> q « 0.6 



22 



GeV is the average energy of the light quark, which coincides with 1/2 of the denominator in 7, while £/[|3 from Table 
IVIIl is iq n &ikn- Thus indeed one has equality M u 7/123 = 7 2/123- 

For practical reasons we have used for our calculations the reduced (2 x 2) forms everywhere. 

TABLE VIII: Dirac eigenvalues £0 (in GeV) for quarks of different masses m (in GeV) and a a . The averaged potentials (U), 
(V) (in GeV) for different a s are also presented. 



Cts 







0.3 


0.39 


m = 


0.005 


0.65 


0.493 


0.424 


m = 


0.15 


0.80 


0.584 


0.509 


m = 


0.2 


0.838 


0.617 


0.539 


(U) 




0.573 


0.486 


0.463 


(V) 







-0.198 


-0.273 



In the nonrelativistic limit one has 



and for the plane- wave (free) quarks, v = e , u Jf ' > , one has v c = — 



C s b — iM u v c ^r-vd i u (C.3) 
J 2m 

■ku„, , . _ - „_;ku 



C sb =i^^. (C.4) 
Am 

Appendix D 
DERIVATION OF EQ.(T40D etc. 

To introduce the Weinberg method it is useful to start from the well-known Hilbert-Schmidt method in integral 
equations with symmetric kernels K(x, y), where x, y belong to the n - dimension space. The eigenvalue equation has 
the form 



(j) n [x) =l n J K(x,y)<f> n (y)dy (D.l) 
The spectral decomposition and the resolvent are 

In 

T(x,y;l) = ^ (D.2) 



and the orthonormality conditions: 



(D.3) 



(j> n K{x, y) 4> m dxdy=\- 8 mn (D.4) 



In the case discussed in section ITVl one arrives to Eqs. (l37H4T)]) . starting from equation 
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and performs symmetrization, using definitions cf> n — \/Hq — E 

K = — 1 V— 1 ln = — ■ (D.6) 

v Ho — E v Hq — E t\ v 

Now Eq. (|D.3J) yields (|42|) . where a v is denned in (|H]l . Eq. (|D.4[) gives (|39|) . Similarly, the Greens function is connected 
to the resolvent 



H -E + V v'Ho-E VHo-E 
^ (l-r?„) 



Now we turn to the i-matrix. One has 

i = T>-t>GV"; H = H a + V (D 
where y = Vm in sector I. One can rewrite (|D.8I) 



= H a -E + V(iJ - £) - J— (flo - £0 

V ^ - 1 

a v (p,E) a v (p',E) 



H -E + J2- 



= H -E-T / ! T^ Z -T / a v a u . (D.9) 

For the latter sum one writes 

= Y J (E n - E) c v n c v n , * n = J] (£7„ - E) „ = (p I (fib - £7) , (D.f 0) 

unn' n 

where was used. Hence finally one gets Eq.(36) 



\ - ^ av{p,E) a v {p',E) 



Appendix E 

ANALYTIC STRUCTURE OF WEINBERG AMPLITUDES AND POLE POSITIONS 

In this Appenix we study the analytic structure of production and scattering amplitudes induced by CC resonances. 
We consider two types of amplitudes, the scattering amplitude in the sector II, e.g. A{DD* — » DD*), and production 
amplitude of the type {QQ) — > (Qq)(Qq), which appears in processes e.g. e + e~ — > DD*,... or B — > K(QQ) — > 
K(DD*). 

In the first case the relevant part of amplitude is given in (|40p . and can be written as 

l-r) u (E) l-rj v {E) 
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In the second case one can start from (|D.7[) for (QQ) Green's function and persuade oneself that neglecting mixing 
of states one returns to the expression (|10p . The production crossection is proportional to the imaginary part of Gqq 
on the cut, starting from the threshold of interest (e.g. DD*) and can be written as 

\M E t - 1 AG aa - y. .aw <-> 

One can easily find, that the latter expression is proportional to ^f v (l) R^~7^p ^i/(2), so that of crucial importance 

is the analytic structure of rg^ ■ 

We consider the case, when only one bare QQ state £"„ is retained, assuming that other states are far off and mixing 
of states, discussed in section [V] is unimportant as compared to the direct influence of the decay channel. In this case 
one can write 



and we write w nn (E) = w(E) 
We write w(E) as 



VAE) = f^g (E.3) 



' (2tt)3 £ - £7(p) 2 y z - u ' 1 ; 



where c = £(p) = + |L, z = 2M (£ - and finally 

/(«) - /(p 2 ) = (J(p)) 2 . (E.5) 

Since f(p 2 ) > for all real p 2 , one has 

- rOG j 

w(0) = -- ^%f(u)<0. (E.6) 

It is convenient to continue f(u) analitically in the region near the real axis 4 and rewrite (|E.4j) as a contour integral 
along the contour C circumjacent the cut in the u-plane 



c f y/udufiu) 
w(z) — ' 



4 J c z — u 

It is clear that the same integral along the contour C with the point z inside C 1 does not have singularities on the 
first sheet, hence one can represent w(z) as follows (difference of two integrals is the residue at the pole u = z) 

W (z) = - Y c^f(z) + F(z) (E.7) 

where F(z) is a nonsingular function which can be Taylor expanded around z — 0. 

In (1E.7|) the argument of z is chosen in a standard way: arg(z) — 0, for z = \z\ + iS, and argz — ir for z < 0. 
We turn now to the analytic structure of Weinberg amplitudes, which using (j4"0|) we write as 

A(E) , -1™ = s (E.8) 

l-n{E) z- z p +mMc^zf{z)-2MF{z) 



4 This is always possible in our Gaussian ansatz for wave functions and subsequent Fourier transform J(p), in more general case one 
continues the absorptive part, as it is used in the dispersion relation technic, via the relation Abs (/(E)) = j*(f z (E) - f n (E)), where 
f l (E) is analytic function denned on the i-th Riemann sheet. In the general case one might encounter potential type singularities in 
complex plane, separated from the positive real axis 
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where we have defined z p = AM (E p — E t h) and E p is the bare position of the QQ level. The denominator in (IE.8[) 
can be rewritten as 

D{z) = z-z p + ib^fz~ /(0) - z p (1 - 77(0)) + n{z) (E.9) 

where we have used relations: 

q(0)= ^(0) = 2Af W (0) ^ 

-E"*?! — Ep £p 

since w(0) = F(0), and 77(0) = - 2M f (0) ■ We also defined b = nMc = Ml, and 

n(z) = «>>/£(/(*) - /(0)) - 2M(F(z) - F(0)). 

Since nO) = 0, we expect it does not affect strongly the analytic structure of D(z) near 2 = 0, where n(z) can be 
written as 

r)(z) = Cl z + i di z 3/2 + 0(z 2 , z 5/2 ). (E.10) 



k = k + , fc_, with 



The poles of A(2£) in the zeroth approximation (n = 0) are easily found, denoting y/z = k, one has two poles at 



+ f (T) a+ ^-^)) 



(E.ll) 



fc_ = +z p (l-» 7 (0)) (E.12) 



i6/(o) / fbf(oy 2 



Here occur two limiting situations, (i) z p is small (the bare pole is in the proximity of the threshold), or 



:„(!-,/«)))•• l^P) (E.13) 



(ii) z p is large, (pole E p far from threshold) 



6/(0) 



2 



|zp(l-7?(0))|» (^) ■ (E.14) 



In case (i) the poles are (neglecting higher order terms) 

~ 6/(0) 



Z„(l-77(0)) 

fc + = -{ PV (E.15) 



One can see, that for weak CC interaction, when 7/(0) < 1, both poles are on the second sheet (virtual states), while 
for strong CC interaction, 77(0) > 1, the pole k + is a bound state, while fc_ is a virtual state. 
Now for the case (ii) one can write 



\ i 6/(0) 



-V(0))J 2 

and in the standard situation, when 21(1 — 77(0)) > 0, one has a pair of Breit-Wigner poles BqT y, with 



(E.17) 



E = El (l- V (0))-^M) I (E.18) 
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r =£eM/(0), Pp = y/2M(E p - E th ) (E.19) 

Note, that (lR19l) coincides with $T5§ as it should be. Using dOf . (fE7T0|) one can write the following analytic 
representation for the Weinberg amplitude in terms of variable k = y^i 

K '> {k-k + )(k-k_)+ Cl k 2 + id 1 k 3 + 0(fc 4 ) V ' ; 

Note, that for the CC poles (k + , fc_ near threshold) the form of A(k) is far from the Breit-Wigner type and is of 
the cusp type, with infinite energy derivative near the pole, which possibly explains the very narrow peak of X(3872). 

Finally, we discuss the case of several thresholds, e.g. in X(3872) for isospin zero one has a sum of DqDq+ h.c. 
and D + D*_ +h.c. terms in w, so that in general case one can write for n thresholds. 



JV 



w[E ) =J2~ T ^^J 1 ^ (E.21) 



i=l 



2 



where Zi = 2Mi {E — E^). One can apply to w(E) the same procedure as before to separate out nonanalytic terms, 
with the result, that D{z) now has the form 

n 

D(z) =z-z p + ib 1 Vzfi(0)+iJ2b i y/z-Aifit-Ai) + n(z) (E.22) 

i=2 

where we have kept notations for z with respect to the lowest threshold, and Aj = 2Mi(E^ — eQ). 

It is important, that for z < Aj the argument of the square root term is (i?) leading to some renormalization of the 
term z p for large Aj, while for small Aj the situation is complicated and should be solved explicitly in the complex 
plane z. 



